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RADAR ALTIMETER AS A SPACEBORNE NAVIGATION AID 

by 

F. Kalli ‘ 

ABSTRACT 

The purpose of this report is to show the feasibility of using a radar 
altimeter as a spaeeborne navigation aid, particularly on maimed space- 
craft by showing, at least in a cursory way: (1) the areas or mission 
phases where it could be advantageous; and (2) that the power, antenna 
size, and other system requirements are both practical and within the 
present state-of-the-art. Although the mission analyses contained herein 
are simplified to permit a cursory evaluation of feasibility, there is suf- 
ficient information to permit an engineer, unversed in orbital mechanics 
to make some "quick -look, desk'* computations. Furthermore, the ap- 
pendices contain some parametric curves for "quick-look” evaluations of 
atmospheric drag effects on near-earth spacecraft. 
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SUMMARY AND CONCLUSIONS 


A radar altimeter could be a valuable navigation aid on-board a manned 
spacecraft, particular whenever the spacecraft is out-of-view of ground support 
facilities, such as 

1. During certain abort trajectories whenever the spacecraft is within the 
useful range of the altimeter (~400 km). 

2. When the spacecraft is at the backside of the moon or a planet. 

3. Possibly, during the non-blackout periods of reentry. 

4. While approaching the moon or planet. 

5. While in orbit about the Earth, Moon or planet. 

For instance, after one orbit about the central body (earth, moon, or planet), 
the spacecraft can use the altitude measurements in the following stepwise 
manner: 
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v„ = / u r /a r 
a ' p a 



where r is the distance from the central body’s center, fi is its gravitational 
constant, R is its radius, h is the spacecraft’s altitude, "a" is the orbit’s semi- 
major axis, e is the orbit’s eccentricity, r and r are the minimum and maxi- 
mum r (in the case of the earth, they are the perigee and apogee radii) and v is 
the spacecraft's speed. 

Assuming altimeter measurement errors of ±100 m (rms, after smoothing), 
normally distributed and uncorrelated, and for an orbit about a spherical central 
body with R = 6367.4 km (the average earth's radius, h p = 200 km, h a = 400 km, 
then the corresponding rms errors in the above determined orbital elements 
would be: 

8 r = 8 r = 8 r = 8 h = 100 m 
p a 

S a = 141 m 
8 e = 1.05 x 10" S 
Sv < 5 8v < Sv ~ 0.12 m/sec 

This leads one to conclude that the spaceborne altimeter could be a powerful navi- 
gation tool, which could also be used to determine the shape of the central body 
if ''a priori" information is available from ground tracking data about the orien- 
tation of the orbit plane relative to the central body's axis. This latter use would 


be highly desirable (see ref. 20) . 
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Regarding the altimeter system requirements, for S/N = 14 db, and effective 
antenna' area of lm 2 , a pulse width of Ip. sec, a pulse repetition rate of 500 pps 
(a duty cycle of 5 x 10 -4 ), an average radar cross-section per unit area of planet 
surface of 0.01, and at an altitude of 300 km, the required peak power output of 
the transmitter is 270 watts, which corresponds to an average power of 135 
milliwatts. Off-the-shelf magnetrons are available with such capabilities. 
Furthermore, by using CW and phase-lock loop techniques to narrow the effective 
bandwidth, either the average power or antenna diameter may be reduced. Thus, 
it is concluded that the radar altimeter could be a practical and useful navigation 
aid, especially for manned spacecraft, and systems trade-off studies are 
warranted. 


x 



RADAR ALTIMETER AS A SPACEBORNE NAVIGATION AID 

I. INTRODUCTION 

The purpose of this report is to show the feasibility of using a radar altimeter 
as a useful navigation aid at altitudes of up to about 200 to 400 km. 

Present day jet aircraft successfully use pulse type radar altimeters at alti- 
tudes of up to 12 km (40,000 ft). It will be shown that by increasing, within 
practical limits, the transmitter power output and antenna gain over and above 
what is presently used on jet aircraft, such an altimeter could be feasibly and 
successfully used at altitudes of up to about 200 to 400 km. Furthermore, the 
following section on mission analysis will show how such an altimeter could be 
used as a valuable navigation aid. 

II. MISSION ANALYSIS 

A radar altimeter could be a valuable navigation aid onboard a manned space- 
craft, particularly during those mission phases of trajectories where the space- 
craft is out-of-view of ground support facilities (stations, ships, or aircraft). 

For instance, it could be useful in the following conditions: 

1. During contingencies such as abort trajectories especially when the 

spacecraft altitude is within the useful range of the altimeter, i.e., ~400 
km. 
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2. During lunar or planetary missions when the spacecraft is on the backside 
of the moon or planet, and the altimeter could be used to give a time 
history of altitude. If the spacecraft gets too close to the lunar or plane- 
tary surface, then corrective actions could be taken such as thrusting 

the spacecraft to a higher altitude. 

3. Perhaps during the non-blackout periods of reentry. Although this area 
will not be discussed further in this report, it merits mentioning because 
of the frequent, and large reentry footprints expected in the future with 
those space laboratory type missions which will utilize lifting body space- 
craft with high lift-to-drag ratios (L/D) of about 2 or greater as ' ’shuttle 
buses" to earth and hack for logistic support. (See refs. 1, 2, 3, 4.) 

4. While in an orbit around the earth, moon, or planets, wherein a time 
history of altitude could be used to determine the three Keplerian orbital 
elements which define the shape of the orbital conic but not its orienta- 
tion, via the following basic equations, wherein first order and higher 
order effects due to gravitational perturbations, atmospheric drag, 

etc., are neglected here for the sake of simplicity. The first order gravi- 
tational effects and atmospheric effects are discussed in Appendices A, B. 

Consider the spacecraft to be in an orbit around the earth, for example. Then 
the spacecraft distance (r) from the earth's center is related to spacecraft 
altitude (h) and the Keplerian orbital elements as follows: 
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( 1 ) 


r(t) =R +h(t) = 


a (1 - e 2 ) 

1 + e cos. 9 (t) 


where as shown in Figure 1 

R e = radius of the earth, strictly speaking at the suh-satellite point 
h (t) = altitude above the earth’s surface at time t 
a = semi-major axis of the Keplerian orbit conic 
e = eccentricity of the Keplerian orbit 
#(t) = true anomaly at time t 


The angle 9 is related to time t in the following, somewhat complicated. 


way: 




(E - e sin E) 


(2) 


where t = time at which the spacecraft passes through perigee 


E = cos -1 



= eccentric anomaly 


(3) 


sin E = 


r . sin 6 


(4) 


tan 


6 _ 

2 


/I +e 

V 1 - e 



(5) 
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This latter equation is perhaps the most useful relation between 8, and E, since 
9/2 and E/2 are always in the same quadrant (see Figure 2.4 /page 39, ref. 5). -At 
minimum altitude (h . ), 9-0 and 

r'(t) . = r = R +h. = a (1 - e) = orbital radius at perigee (6) 

^ / mm p e 1 mm v r v 1 

At maximum altitude h , 6 = tt radians and 

max 

r ( t ) max = r a = R e + = a (1 +e) = orbital radius at .apogee (7) 

The spacecraft speed relative- to the earth's center-of-mass can be determined 
from 



( 8 ) 


, Spacecraf t’s orbital- (9) 

speed at perigee 


Spacecraf t’s orbital (10) 

speed at apogee 

jj. = earth's gravitational constant 
= 3.986 032 x 10 14 m 3 /sec 2 
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Since the earth is not spherical in shape,, and there are mountainous regions, 
a good approximation is obtained by replacing R g with the average radius of the 
earth <(R e )> in the equations above. Using the constants given in reference 6, we 
have 


<X> = 


-7J-/2 

Jo 


R eq (1 - f sin 2 6) dQ 


TT 

~2 


where here Q = geocentric latitude 

R = earth's equatorial radius 
= 6,378,165 km 
= 3443.934 n. mi. 
f = earth's flattening 


( 11 ) 




f = 


R eq — R pole _ 1 


R 


eq 


298.30 


( 12 ) 


R P oie = eart ^' s polar radius 
= 6,356.77 km 


= 3432.38 n. mi. 


Integrating equation 7, 


< R e> = R „ 



= 6367.47 km = 3438.15 n.mi. 


"-(13) 
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Thus, it can be seen from the above discussion that adequate equations exist 


so that a small onboard computer could be used to determine the three Keplerian 
elements a, e, and t . It should be noted that these three elements define the 

P 

conic but not the orientation of the conics plane in space, wherein it is assumed 
that there is no ’’a-priori 1 ' knowledge of the spacecraft's position vector or 
velocity. 

Let us now digress momentarily and examine roughly how the altimeter’s 
accuracy affects its capabilities as a spaceborne navigation aid. Performing a 
regression of h on t to get a least squares curve fit for the data and to thence 
get a "band" or "error" for h min is beyond the scope of this report. Further- 
more, such a sophisticated approach is not warranted at this time because of 
the many simplifying assumptions that are used. Hence, the following simpli- 
fied analyses are only cursory, and their only purpose is to show feasibility. 
Assuming for simplicity that the spacecraft has completed one orbital revolu- 
tion about a spherical body of radius It, then 


r (t) = h (t) + R 

(14) 

r = h . +R 

p mm 

(15) 

r = h + R 

a max 

(16) 

r + r 

a p 

a = p 

2 

(17) 

r — r 

e- a p . 

(18) 
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( 19 ) 


V 

p 



V 


a 



( 20 ) 


v(t) - 



( 21 ) 


Neglecting variations in the surface terrain of the central body for simplicity, 
then the errors in r and r are 

P - a 


Sr = Sh . = Sh) the altimeter’s read-out accuracy after smoothing 

p min 1 J a 


Sr = Sh „ = Sh 

a max 


Sr = Sh 


( 22 ) 

( 23 ) 

( 24 ) 


Assuming a normal distribution in the altimeter measurement errors and they 
are uncorrelated 


Sa = 



( 25 ) 


Sa = /2 Sh 


( 26 ) 
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In a similar manner, we obtain for the error in e , 



Se % Sh, for 0 5 e ^ 0.2 

2 a 


while for the error in v, we get 


( 27 ) 


(28) 


(29) 

— = v for 0 5 e ^ 0.2 (30) 

v r 2 r 

These latter equations are applicable for any point in the orbit, i.e., they can be 
used, to determine Sv p or Sv a by merely substituting r p or r a for r, respectively. 

Consider the case where Sh = 100 m (an assumed rms value after smoothing), 
r p = 200 km + R, r a = 400'km + R, and assume that R = 6367.47 dan, which is the 
average earth's radius. Then 

r = 6567.47 km 

p 

r = 6767.47 km 

a 
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a = 6667.47 km 


e = 0.0X5 

Sr = Sr = Sr = Sh = 100 m (rms) 

pa a / 

Sa = 141 m (rms) 

Se = 1.05 x 10" 5 (rms) 

Since Sv p represents the worst case, i.e., Sv is maximum at perigee for a 
given Sh, we will look at Sv p : 


Sv = 0.12 m ■ (rms) 

Thus , it can be seen from the above example, that the radar altimeter could be 
a very powerful tool as a navigation aid. Furthermore, if the spacecraft had 
some "a-priori" knowledge about the orientation of its orbital plane relative to 
the planet's axis, then the altitude data could be useful in determining the shape 
of the planet. 

The reader is cautioned that the above analyses were simplified to demon- 
strate feasibility and hence are only cursory. For instance a spherical central 
body was assumed, whereas the earth is an oblate spheroid to first order and 
not a smooth spheroid at that because of its mountain ranges, etc. Furthermore, 
atmospheric drag effects were neglected. Appendices A and B show to first 
order how the earth's oblate shape (neglecting the effect of mountain ranges) and- 
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atmospheric drag, respectively, affect the spacecrafts altitude. Thus, it can 
readily be seen that further analyses are warranted. However, it is believed that 
such detailed analyses would be so cumbersome that they would detract from 
the intended purpose of this report, and hence should be the subject of a separate 
report. 


m. SYSTEM ANALYSIS 

Let us now examine the altimeter from the viewpoint of transmitter power 
and antenna gain requirements. 

The two-way range equation fori a pulse type radio altimeter is derived as 
follows; 


The incremental signal power (dS gnd ) received by the increment of ground 
area dA gnd (see Figure 2) is 


dS 


P t G t F t (6>, 4 >) | d~A | 


gnd 


4tt- r 2 L, 


(31) 


where P = transmitter output power 


4 Tt Aj, ef f 4 7T 

G = — =7^ 

\ 2 e 2 


(32) 


G t = transmitter antenna power gain 
A„ ,, — effective transmitter antenna area 

T, e f i 
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F t (6, 0) = transmitting antenna radiation pattern 
\ = carrier wavelength 
•jju = antenna radiating efficiency 
= 0.5 for a circular aperture 
0 T ~ transmitted beamwidth at the half power points 
L T = transmitter line losses 

I dA | = component of dA gnd in the direction of incident radiation 
| <dA | = dA gnd cos 6 


h 2 tan 6 d 6 d 
cos Q 


(33) 


h = altitude of spacecraft altimeter above earth's surface 
r = distance of dA . from altimeter antenna 

gnd 


r 


h 

cos 6 


(34) 


Let ct 0 be the radar cross-section (or scatter cross-section) per unit area 
of illuminated ground. By definition the radar cross-section (or scatter cross- 
section) of a target is the area intercepting! that amount of power which, when 
scattered equally in all directions, produces an echo at the radar equal to that 
from the target. In other words 


u - 4 77 P° wer scattered back at the source/unit solid angle 
Power/unit area striking target 


(35) 


In general, cr 0 is a function of 9 , so that henceforth it will be written as o- Q (8). 
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The signal power received by the altimeter at the receiver input terminals 
from the target dA gnd is 


dS = 


dS 


gnd 

47 rr 2 L„ 


-o(^) 


( 36 ) 


where L_ is the altimeter receiver line loss, and A„(d , d>) is the effective receiv- 

R In ' 

ing aperture perpendicular to r and is given by 


A R (0,dO 


G r F r (^ i 0)X 2 

Att 


( 37 ) 


where G R is the receiving antenna power gain, and F R {6 , 4 > ) is its radiation 
pattern. 

Assuming that the altimeter transmitting and receiving apertures are the 
same, then 


dS - 


P T G 2 F 2 <p) A 2 <r 0 i e ) cos 6 dA 


gnd 


64 t r 3 r 4 L T Lg 


( 38 ) 


and the total signal power is 


S = 


P T G 2 X 2 


64 tt 3 L T Lg J A 


l 


F 2 (6, 4>) cr Q (O') cos 6 dA gnd 


( 39 ) 


gnd 
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but, as can be seen from Figure 2, 


h 

cos 9 


(40) 


dA 


gnd 


__ h 2 tan 8 dddcp 
cos 2 6 


(41) 


so that 


S = 


P T G 2 


64 7r 3 h 2 L t L e 


U. 


F 2 ( 6 , 4>) cr Q ( 6 ) sin 6 cos 6 d 6 d<£ (42) 


Shown, in Figures 3, 4, and 5 (see refs. 7, 8, 9, 10) are some radar cross- 
sections for roads, forrest and sea at X-band (10 GHz) as a function of y , the 
angle of incidence (sometimes called grazing incidence or grazing angle), where 
y is the complement of the angle 6/2 used in the above equations and shown in 
Figure 2, i.e,, y = 90 - 9/2. Figure 5 also shows the frequency dependence 
of <r 0 . Thus, it can be seen that the evaluation of the above integral can be 

f 

difficult except in some simple cases . The above expression is valid for C-W 
or pulse type altimeters which are beamwidth limited, that is the beamwidth 
is narrow enough that the difference in the two-way time delay between a ray 
along the center of the beam and a ray along the edge (or half power point) of the 
beam is less than the pulsewidth (r). Thus the radius (■i) of the maximum effective 
ground area illuminated by the beam is from Pythogorean’s theorem (see Figure 6) 
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s l/ hCT + (t ) 2 


Since 


But 


— « her 
2 

■l ^ /her 



and 



( 44 ) 


(45) 


(46) 


where is the beamwidth at the half power points of the radiated beam, and 
is the directional gain and is related to the power gain G by 




= radiating efficiency 


h 


% CrG -vCrG 


mast 


7 T 


(47) 


(48) 


where here h max is the maximum altitude for beamwidth limited operation for a 
pulse type altimeter . 
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Let us now evaluate the equation 21 for a simple case only -to get a feeling 
for the magnitude of the so-called loop loss (S/P T ). Assume the beamwidth is 
narrow with a constant intensity, i.e., assume beamwidth limited operation, 
cr 0 (<9) is constant over the narrow beamwidth, and F (0 ; 4 >) = 1. The loop loss 
becomes 

s _ G 2 \ 2 cr 2 0 
P t 32 7T 2 h 2 

where /2 is the half beamwidth for the case being considered. 

By definition the directive gain is 


i 2 

1 - cos — 
2 


(49) 


477 P 


G D = 


max 


I 


p (0, <£) aa 


or in words 


(50) 


r - maximum radiation intensity 
average radiation intensity 

q _ 47 r (Maximum power radiated/unit solid angle) 
D total power radiated 

but 

P(d,0)=F(d,0)P max 
dfl = sin 6 d 6 d 


(51) 
i (52) 
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and F (9 ,4>) = 1 for the case considered. Integrating equation 50, we .have 


and 



2 


1 - COS 



S t) r G X 2, Op 
16 7 T 2 h 2 


1 + COS 


ft 


B 


For small beamwidths 


( 53 ) 


(54) 


a 


cos 


B ~ 


1 + COS 






(55) 


Using 


G = 


4tt A 


T.eff 


X 2 


(56) 


S ^ ^T. eff °0 

477h 2 


(57) 


This formula is accurate enough for the purposes of this paper and agrees with 
the results obtained in references 10, 11. This equation is valid for the C-W or 
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beamwidth limited pulse type radar altimeters. It is interesting to note that, 
except for the dependence of on frequency, the so-called loop loss is inde- 
pendent of frequency as long as beamwidth limited operation is maintained . 

Let us now examine the transmitter power requirements for a particular 
ease. Consider the altimeter to be used on a spacecraft as a navigation aid for 
determining the energy dependent orbital elements, i.e., the in-plane or 
Keplerian elements, as discussed earlier in the paper. Assume or consider 
the following: 

T = 1 / j . sec 

B = = 10 6 Hz, effective receiver beamwidth (58) 


— = 30 = 14.8 db, ratio of signal power to noise power 


f = 10 GHz, carrier frequency 
<x 0 = 0.01 = -20 db 
N r = 0.5 = -3 db 

L t = L r = 2 = 3 db, to be conservative 
h = 300 km 

Now the noise power N is 


(59) 
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where 


K = 1.38 x 10" 23 watt-sec/deg - Boltzmann's constant, 

T f f is the effective receiver noise temperature at the receiver input terminals 
and is given by 


T ,f, = T„ (F„ - 1) + 



(60) 


where T 0 = 290 °K, reference temperature 

F 0 = receiver noise figure, and for a receiver with a parametric ampli- 
fier, a typical value is about 4 db = 2.5 (see ref. 12) . 

The T a is the antenna noise temperature, sometimes called the sky noise tem- 
perature, and is due to the background radiation, side lobe noise, spillover, and 
atmospheric attenuation. For a receiver looking at the earth with a beamwidth 
equal to or less than the earth’s diameter, T a = 300°K (see ref. 13). The T L is 
the actual temperature of the receiver input terminals and will be assumed to 
be 290°K for a spaceborne receiver. 

Thus, T eff reduces to 

T ff -T 0 F 0 (61) 

£ 725°K 


and the required peak transmitter power is from equation 37 


18 



(S/N) 47rh 2 KT ff BL t L r 


P t = 


^ A 


T.eff ^0 


30 (4ir) (300 km) 2 (l.38 x 10" 23 — .- £ e £ . ) (725 °K) (10 6 Hz) (2) (2) 

\ ± L (62) 

0.5 (1 m 2 ) (10“ 6 Km 2 /m 2 ) (10“ 2 ) 


= 271 watts, peak power 


See also Table 1 which gives the above analysis in logarithmic form. 

Table 1 

Analysis of Altimeter Transmitter Peak Power Requirements for h = 300 km 


Signal-to-Noise Power Ratio 

' S/N = 30 

14.8 db 

l/Spreading loss 

47rh 2 = 1.13 x 10 6 (km) 2 

60.5 db (km) 2 

System Noise Temperature 

KT ,, = 10" 20 W 

ex f 

-200.0 dbW/Hz 

• Noise bandwidth 

B = 10 6 Hz 

60 dbHz 

Altimeter transmission loss 

L T = 2 

3 db 

Altimeter reception loss 

II 

to 

3 db 

l/(radiating efficiency) 

JL= 2 

N r 

3 db 

Ineffective aperture) 

-1— = 10 6 (km)- 2 

^eff 

60 db/km 2 

l/(ground scatter cross- 

JL= 10 2 

20 db 

section per unit area) 


- 

Transmitter Peak Power* 
required 

P T = 271 W 

24.3 dbW 


*At a duty cycle of 5 X 10~ 4 , this corresponds to an average power of 135 milliwatts (see following 
analysis and Table 2). 
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To check if this example is beamwidth limited, we now compute the maxi- 
mum altitude for beamwidth limited operation from 


h 


max 


~CtG 


where 


G = 


47T A 


eff 


A . 2 



3 x 10 8 m/sec 
10 10 cps 


= 3 x 10" 2 m 


G = 477 1 1 ^ _ = 1.4 x 10 4 = 41.4 db 

(3 x iO- 2 m) 2 • 


h 


max 


3 x 10 s km/ sec (IQ- 6 sec) 1.4 x 10 4 
4 


= 1.05 x 10 3 km 

Hence for operational altitudes below 1,000 km this system is beamwidth limited. 

Let us now examine roughly the pointing requirements, because they would 
affect the attitude control requirements. Since 

A ^ o 5 n d2 

A eff - u :> — 4 — 


d = 


8 A eff ' 
77 


1/2 


= 1.6m 
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then 


ft, =— = 3 x 10 2m -= 1.87 x 10' 2 rad 
B d 1.6 ra 

^ B = 1-7 deg 

This means that the altimeter must be pointed in the direction of the vertical axis 
towards the ground (i.e. in the direction of the negative zenith axis) with an 
accuracy of about ±<9 Q /2 or ±0.85 deg. which is well within the state-of-the-art 
even for manned spacecraft. Otherwise the systems gain would be significantly 
reduced whenever the return signal with the minimum time delay, i.e. the signal 
along the altitude axis, would be outside the half power beamwidth because of 
improper antenna orientation. 

Perhaps it should also be pointed out here that as the altitude varies, then the 
Doppler information could also be utilized for navigation purposes, i.e. the 
two-way Doppler shift Af = 2h/k , where h is the time rate of altitude change, 
and X is the wavelength. 

Regarding the problems which may be encountered with uneven terrain, 
particularly on close approaches, both a wide beam and narrow beam could be 
advantageously used. The wide beam would provide information regarding the 
average terrain while the narrow beam would correlate the target being approached 
with the average terrain. 

It is interesting to note that the Apollo, lunar module, landing radar 'uses 
the principles discussed, i.e. it has one altitude sensor beam used also as a 
velocity sensor and three velocity sensor beams. Also, the Surveyor landing 
radar used the above principles. 
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The maximum pulse repetition rate for resolving altitude (range) ambiguities 
without having to resort to coding the pulses is derived as follows: 

At altitude (h) of 300 km, the two-way propagation time (2 t prop ) 

2 t - = 2 x 10“ 3 sec (63) 

prop n 

* 

To avoid the necessity of coding the pulses in order to resolve altitude ambiguities, 
this two-way'propagation time must be equal to or less than the pulse repetition 
period (T r ), i.e. 

2 1 = 2 x 10- 3 sec < T (64) 

prop r v 1 

The pulse repetition rate (f ) is* 


f f = J_= 500 pps. (65) . 

and the duty cycle is r f = 5 x 10“ 4 . 

Shown in Table 2 are some beacon magnetrons which can more than ade- 
quately meet the peak power requirement of 270 W and the duty cycle of 5 x 10“ 4 , 
which corresponds to an average power of only 135 milliwatts. Thus, it may 
be concluded that using an altimeter onboard a spacecraft as a navigation aid - 
is feasible. 
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Table 2 


Beacon Magnetrons 


FIXED FREQUENCY 

All units positive-pulsed except those indicated by an asterisk. 

Model 
Number 
(ref. 14) 

Freq. 

GHz 

Min. 

Peak 

Power 

w 

Duty 

Cycle 

Pulse 

Width 

/xs 

Peak 

Anode 

Volt. 

kV 

Peak 

Anode 

Cur. 

A 

Approx. 

Weight 

oz. 

MA-259 

7.0-7.5 


0.005 


2.00 

mmm 

H&EH 

21 

MA-221B 

8.0-8. 8 


CW 

- 

0.45 

mmm 

11 

MA-231B 

8. 0-8. 8 

20 

CW 

— 

0.93 

0.06 

11 

MA-221A 

8. 0-8. 8 

10 


1.0 

|E1 

0.15 

11 

MA-221C 

8. 0-8. 8 

10 

1 

5.0 


0.15 

11 

MA-221D 

8. 0-8. 8 

20 

1 

1,0 

B 

0.30 

11 

MA-231A 

8. 0-8. 8 

200 

1 

■IB 

1.00 

WSM 

11 

MA-250 

8. 5-9. 6 , 

4kw 


1 

3.75 

isa 

21 

MA-212B 

8. 8-9.6 

1.5 

CW 


0.46 

B 

11 

.MA-212A 

8. 8-9.6 

— 

!0.01 

i.o ; 

0.52 

B 

11 

MA-232B 

8. 8-9.6 

mm 

- CW 

I _ 

0.93 

■SB 

11 

MA-212C 

8. 8-9.6 

B 

0.01 

1 5.0 1 

iQ.52 ■ 

0.15 

11 

MA-212D 

8. 8-9.6 

— 

1 


0.56 

B 

11 

MA-232A 

8. 8-9.6 



BSJ11 

1.00 


11 

MA-252A 

8. 8-9.6 




1.15 

B 

11 

MA-222* 

9. 3-9.4 

1 

0.002 


5.50 

4.50 

48 

MA-239 

13.28-13.37 

mm 

0.25 


0.85 

0.24 

23 

MA-239B , 

13.0-13.6 

B 

CW 

B 

0.95 

0.06 

23 

MA-240 

13.7-14.6 

500 


mm 



16 

MA-240A 

13.7-14.6 

700 


m 

1.75 

1 

16 

MA-245 

15.4-15.6 



B 



23 

MA-257 

16.0-16.5 

lkw 

0.002 

0.5 

2.50 

'1.80 

16 

MA-256* 

16.4-16,6 

lOkw 

0.001 

0,2 

9.00 

5.25 ( 

22- 

MA-246A 

16.5-16.8 

100 

0.005 

0.5 

1.15 

0.45 

16 


Negative pulsed 
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Table 2 (Continued) 


TUNABLE 

All units positive-pulsed except those indicated by an asterisk. 



Min. 


Pulse 

Width 

jJL S 

Peak 

Peak 

Approx. 

Weight 

oz. 

Model 

Number 

Freq. 

GHz 

Peak 

Power 

w 

Duty 

Cycle 

Anode 

Volt. 

kV 

Anode 

Cur. 

A 

MA-219F 

8.5-9.0 

3 

CW 

_ 

■as 

wtrm 

13 . 

MA-219B 

8.5-9. 6 

1 

CW 

— 

BSE 

| 

13 

MA-214H 

9.0-9. 5 

3 

CW 

— 

0.65 

0.02 

13 

MA-214B 

9.0-10.0 

1 


— 

0.47 

0.02. 

13 

MA-249 

9.25-9.45 

17 • 



ns 

0.05 

13 

MA-219A 

8.5-9.6 . 

10 




1.0 

mm 

0.15 

13 

MA-219C 

8. 5-9. 6 




0.56 

Ed 

13 

MA-261 

8. 5-9. 6 

50 



0.82 

MEM 

13 


8. 9-9.4 

400 

0.0005 


4.00 

0.5 

24 


8. 9-9.4 

900 



4.30 

0.8 

24 


8. 9-9.4 

800 


I 

4.30 

0.9 

24 

MA-214C 

9.0-10.0 

8 



0.54 

0.15 

13 

MA-214D 

9.0-10.0 

16 


1.0 

0.57 


13 

MA-241 

9.1-9.4 


n 

1.0 

2.75 


24 

MA-218* 

9,3-10.0 

7kw 


1.0 

5.90 ' 


48 

MA-232T 

9. 5-9. 8 


0.002 

1.0 

Id 

0.7 

13 

MA-260 

16.0-16.5 

lkw 

0.001 

0.5 

B9 

1.6 

18 


(per ref. 14) 


If C-W operation is desired because a transistorized transmitter can be used 
without the need for filament power, then the average power output required is 


P ave = p T Tf r = 270 (lO' 6 ) (500) = 135 mW 


for the example being considered. 














































However, the bandwidth must also be reduced by an equivalent amount in order 
for the sample system being considered to operate at 300 km altitude and S/N = 
= 14.8 db. That is 


B cw = Bpuise T f , = 500 Hz (67) 

Using phaseTlock loop techniques, this bandwidth is not only feasible but can be 
reduced an order of magnitude, in which case the above corresponding average 
power requirement can also be reduced an order] of magnitude. However, it may 
be more desirable to reduce the antenna aperture instead of the average power. 

In any event, system trade-offs are possible, and the altimeter is a feasible navi- 
gation aid for use on spacecraft. Furthermore, as earlier analyses in this paper 
indicates, the altimeter's measurements could be useful in determining the 

shape of the moon or planet, if "a-priori" knowledge is available from ground 
I 

.tracking data about the orientation of the spacecraft's orbital plane relative to 
the axis of the central body. Determining the shape of the moon and planets is 
one of the significant goals of the scientific community (see ref. 20). 
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APPENDIX A 


Earth Gravitational Perturbations to First Order 

A near-earth satellite experiences altitude variations because of gravitational 
perturbations due to the non-spherical shape of the earth; i.e., the satellite orbit 
is not a true ellipse. The satellite orbit to first order is given by (ref. 15) 


r 0 1 - e cos (J3 - cS) — — 

6 V r o 


sin 2 i 0 cos 2 fi 


whenever the eccentricity e is of the order of J , where F Q is a constant of 
integration and has the physical significance of being essentially the mean, semi- 
major axis of a perturbed orbit, J = 1.6235 x 10“ 3 is the coupling coefficient 
of the second harmonic in the earth’s gravitational potential (J = 3/2 j 2 ), R eq = 
6, 378.165 km (3,433.934 n. mi.) is the earth’s equatorial radius, i 0 is the mean 
orbital inclination about which the instantaneous inclination varies, o> is the 
angle of perigee, and J3 is the central angle. (See Figure A-l which depicts the 
orbital geometry.) 

The satellite altitude h is 

h = r - R e (A-2) 

where the earth’s local radius R e is given by, to first order 

R e = R eq (1 - f sin 2 d) 

. 29 


(A-3) 



and 


f = — — = earth’s flattening 

298.30 

9 - geocentric latitude 

sin 6 - sin i Q sin /S (A-4) 

For the case of e = 0, F Q = R eq + 185.2 km (i.e. h = 100 n. mi. at the equator), 
and i 0 = 90° (a polar orbit where the altitude variation are maximum; see 
equation A-l), the altitude varies sinusoidally about an average (h) of ~ 200 km 
with an amplitude of 0.8 km as shown in Figure A-2. It is interesting to note 
that the earth's radius varies by 21.50 km from equator to pole, see Figure A-3. 
Hence, it may be concluded that in this example, the’ satellite is closely following 
the oblate shape of the earth. 
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APPENDIX B 


Atmospheric Drag 

A near earth satellite experiences a dissipative drag force as it goes through 
the earth’s atmosphere. The drag force per unit satellite mass is 

F - -B/jvv (B-l) 

where B is the satellites ballistic coefficient defined by 


2m 


(B-2) 


Cp is the drag coefficient which depends on the satellite shape, its orientation as 
it moves through the air, the molecular density of the atmosphere, and the 
satellite's air speed (v). (See ref. 16). For massive, near-earth satellites such 

# I 

as Apollo, for instance, the = 2.0 ± 0.2 (see ref. 6). The A is the satellite 
cross-sectional area normal to the satellite air velocity v , and m is the satellite 
mass. 

This drag force is called a dissipative force because it reduces the total 
energy of the satellite via the air friction. To illustrate this, consider first for 
simplicity the case of a near-earth satellite in a circular orbit. Assume for this 
illustration that the earth’s atmosphere is spherical in shape and is not rotating 
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(the more general case is discussed below later on). The satellite's kinetic 
energy (T) per unit mass is v 2 /2, and its potential energy (U) per unit mass is 
-p / r , where p is the gravitational constant of' the earth (p- 3.986 032 x io 14 
m 3 /sec 2 ; ref. 6). The satellite's total energy (E) per unit mass is 

E = T + U (B3) 


In a Keplerian orbit 



(B4) 


•where 

a = semimajor axis of the elliptic orbit 

Substituting from Equation B4 into Equation B3, the satellite's total energy 
per unit mass becomes 

E = - — , elliptic orbit (B5) 

2 a 

In the case of a circular orbit (a = r ) , 

circular orbit (B6) 

2 r 

The atmospheric drag changes the satellite's total energy by an amount W D , the 
work done by the drag force, i.e. 
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The amount of change in total energy in one orbital revolution is 


AE 

rev 



(B8) 


and also, from Equation B6 


AE _ [x Ar 
rev 2 r 2 rev 


(B9) 


Since F D is antiparallel to dr, and is essentially constant over one orbital revolu- 
tion for a circular orbit in a spherical atmosphere 


I 


dr = K 


1 

‘'rev 


dr = - 2 7rr B p vl 


(BIO) 


where v c is the circular orbit speed and is given by 


v 


2 _ p- 

c — 


(Bll) 


so that 


AE 


2rrB p a 


rev 


(B12) 



The minus sign here shows that the dissipative drag force decreases the total 
energy of the satellite. Combining Equations B9 and B12, and solving for Ah/rev 
which is equal to Ar/ rev. 


Ah 

rev 


4-77-Bpr 2 , for circular orbit, and a circular 
non-rotating atmosphere 


(B13) 


The minus sign merely indicates that the altitude is decreasing due to atmospheric 
drag. 

The case for an elliptic orbit in an oblate rotating atmosphere has been 
treated for the cases of 0 5 e 5 0.01 (ref. 17), and e - 0.01 (ref. 18), and will 
not be treated here except to make the following observations. It is interesting 
to note that both the rotation of the atmosphere and its oblate spheroidal shape 
contribute significantly to the orbital decay. For instance, for a circular orbit 
in the earth* s oblate, rotating atmosphere 

— = - 4tt B r 2 (1 _ d) 2 I (q/2) P n e' q/2 (B14) 

rev 

\ 

where I 0 (q/2) is the Modified Bessel function of the first kind; I Q (q/2) e" q/2 
accounts for the atmosphere's oblateness; (1 -d) accounts for the atmosphere's 
rotation; 


p n = atmosphere density at the minimum altitude (h min ) of the. satellite, 
q = KR eq f sin 2 i Q (B15) 
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(B16) 


p ~ e Kz , atmosphere density at altitude h 
K = - (d/dz) -£n p, lapse rate of the atmosphere (see ref. 19) 

f 

e = base of the natural logarithms 

z = height above the minimum altitude (h . ) 

f = 1/298.30, earths flattening 

R eq = 6,378.165 km, earth’s equatorial radius 

d =Q S /n (1 - e 2 ) 1/2 cos i 0 (B17) 

e = orbital eccentricity; and is zero in the case for the circular orbit 
being considered. 

= 4.178 07416x10 _ 3 deg/sec, -earth's rotational speed about its axis 
n = 2tt/t, satellite's mean motion (B18) 

r = anomalistic period, i.e. the time in orbit from perigee to perigee 
i 0 - orbital inclination 
For near-earth orbits and e < 0.01, 

d =^-cos i 0 (B19) 

From the above equations, it can be seen that: the effect of the atmosphere’s 
rotation on the satellite's altitude decay (Al/rev), due to the (1 - d ) 2 term, 
could be as large as ± 13% per orbital revolution and is greatest when cos i Q = 1; 
while the effect of the atmosphere’s oblateness on the satellite's altitude decay ' 
(Ah/rev) - , due to the I Q (q/2) £~ q/2 term, could be as large as -25% at satellite 
altitudes of 200 km and is greatest when sin i = 1. Hence, both the oblate shape 
and the rotation of the atmosphere are significant contributors to the satellite's 


altitude decay. 
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(B20) 


(B21) 


(B22) 


For circular earth orbits and a nominal atmosphere, Figure B1 shows the 
ratio Ah/ rev/B ftf f as a function of satellite altitude, while Figure B2 shows the 
ratio B gf f /B actual as a function of both altitude and inclination angle. To use 
these figures, one must know the satellites actual ballistic coefficient, i.e. 

C^A/ 2m, the circular orbital altitude, and the orbital inclination. Then one may 
obtain: 

(a) the B eff from the B eff /B actoal given in Figure B2. 

(b) the Ah/rev from the Ah /rev/B ef f given in Figure Bl. 

It should be pointed out that in the above analysis and considerations, a 
Keplerian type orbit was assumed. What needs to be considered further ,is both 
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the gravitational perturbations (at least to first order) and the oblate, rotating 
atmosphere simultaneously acting on the satellite's motion in order to derive 
a closed form solution of the type given above. 
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SPACECRAFT AT TIME t 



88 



ZENITH 

I 



dA 


= dAgnd cosd = 


h 2 tan# 66 dft 
cos0 


Figure 2-Geometry for Computing Effective Ground Reradiating Area. 
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Figure 3-Radar Cross-section (o^) of Gravel and Concrete Roads at X-Band 
(10 GHz) versus Angle of Incidence ('>'). 
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AVERAGE FORREST CROSS SECTION 
PER UNIT AREA (a°>, db 



Figure 4— Average Radar Cross-section per Unit Area of Forrest (Oq) as a 
Function of Incidence Angle (y) at X-Band (10 GHz). 
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AVERAGE RADAR CROSS SECTION 
PER UNIT AREA OF SEA (a°), db 



ANGLE OF INCIDENCE OR GRAZING ANGLE (y), DEG. 

Figure 5— Average Radar Cross-section per Unit Area of Sea as a Function of 
Wavelength and Incidence Angle ( 7 ). 
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ALTIMETER ANTENNA 



h = ALTITUDE 

9 b =beamwidth at half power points 


Figure 6— Geometry for Computing the Maximum Altitude for Beamwidth Limited 
Operation of a Pulsed Type Altimeter, 
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EARTH SATELLITE ALTITUDE, km 



CENTRAL ANGLE (g), DEG. 

Figure A2-Satellite Altitude Variations About an Oblate Earth to First Order; 

e = 0, i = 90°, h'= 200 km. 
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Q I ! I 2 — I ^ I L_ 

— 90 -60 -30 0 '30 60 


EARTH'S GEOCENTRIC LATITUDE (0), DEG. 

Figure A3— Earth's Radius as a Function of Latitude to First Order. 
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EVOLUTION PER EFI 





RATIO OF EFFECTIVE TO 
ACTUAL BALLISTIC COEFFICIENT 



INCLINATION ANGLE, DEG. 

Figure 82-Ratio of Effective to Actual Ballistic Coefficient for Circular Orbits 
in the Nominal Earth’s Oblate Rotating Atmosphere. 




